Specific heat of underdoped cuprates: RVB versus Fermi arcs 
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A recent microscopic model of the pseudogap state, based on the resonating valence bond (RVB) 
spin liquid, has provided a simple ansatz for the electronic self energy in which a gap forms on the 
antiferromagnetic Brillouin zone as the limit of a Mott insulator is approached in the underdoped 
regime. Here, the ansatz is employed to calculate the electronic specific heat when a superconduct- 
ing gap is also included. We find qualitative agreement with all experimental observations in the 
underdoped regime of the cuprates. We explore the relationship of the theory to two other purely 
phenomenological approaches, the nodal liquid and the Fermi arc model, and provide justification 
for their use on experimental data in light of this microscopic RVB theory. 
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Simplified models, such as BCS theory, which ignore 
much of the complicated details associated with actual 
metals, have had an enormous impact on our understand- 
ing of the superconducting state. A generalization of 
BCS to include the d-wave symmetry of the supercon- 
ducting gap in the cuprates has taken us a long way in 
understanding the overdoped and optimally doped part 
of their phase diagram. However, it has failed so far to 
provide even a first qualitative picture of the observed 
properties as the Mott transition is approached in the 
underdoped regime. Some additional essential element 
is still missing which has been widely recognized as as- 
sociated with the opening of a pseudogap, although the 
exact nature of the phenomenon remains unknown. Re- 
cently, Yang, Rice and Zhang^ have provided a simple 
model for the self energy in the pseudogap phase from 
which the coherent part of the electronic Green's func- 
tion can be constructed. The work is based on results for 
a resonating valence bond spin liquid and contains a gap 
which is formed on the antiferromagnetic Brillouin zone 
(AFBZ) as the doping is reduced and the transition to a 
correlation-induced Mott insulating state is approached 
from the metallic side. Since its appearance in 2006, 
the model has had considerable success in understand- 
ing some aspects of Raman- and optical properties^ of 
the underdoped cuprates and has also been applied to 
angle-resolved photoemission (ARPES) data,"* 

In view of these developments, it is very important 
to test the YRZ model on other data. In this paper, 
we consider the electronic specific heat, which has long 
been known to show anomalous properties^i^ii not de- 
scribable within simple BCS theory. Lacking the ex- 
istence of a well-developed and accepted microscopic 
theoretical framework, Loram et ali^^^il included a de- 
pression of the electronic density of states (DOS) near 
the Fermi energy to analyze their early specific heat 
work. More recently, a similar analysis applied to optical 
data®'^ has also yielded new insights, including clarifica- 
tion of temperature-dependent Fermi arcs. More sophis- 
ticated, but still purely phenomenological, approaches to 
the pseudogap phase have appeared, including the idea 
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FIG. 1: (Color online) (a) Plot of 7/70 vs T/to for the YRZ 
model with x = 0.16, 0.12 and 0.09. (b) Plot of normal- 
ized jump, AC/Cn ~ AC /'y{Tc)Tc, and normalized internal 
energy, AU/AUn as functions of doping, x. 



of a nodal liquid^^ and of temperature-dependent Fermi 
g^j^^ gio, 11^12, 13^14 ^ In these models, the pseudogap exhibits 
d-wave symmetry and forms on the Fermi surface. At 
high temperatures only the antinodal region is gapped 
while the remaining ungapped arc length about the nodal 
direction is proportional to temperature. The limit of 
a nodal liquid is when the pseudogap is taken to form 
over the entire Fermi surfacej^Siil In comparison to these 
models, the YRZ approach is more securely based in mi- 
croscopic theory. It is also profoundly different in that 
the pseudogap forms on the AFBZ. After presenting our 
results, we will provide an analysis of how YRZ relates 
to both the nodal liquid and Fermi arc model. 

There are of course many other theoretical approaches 
to the pseudogap phase, for example, preformed paira^^, 
existing below an energy scale set by the onset temper- 
ature, T*, with phase coherence taking hold only at the 
lower superconducting Tc, or as extensions of BCS theory 
to include the formation of finite momentum pairs which 
persist above TcJ^ Another class of theories involves 
competing orders, such as d-density wavesii, which set in 
at T* and can coexist with superconducting order below 

In the YRZ model, both the superconducting gap. Age, 
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and the pseudogap, Apg, have a d-wave fc-space depen- 

dence described by: Age — -^(cosfc^a — cos/cj,a) and 

Apg = — ^(cosfcj.a — cosfcj,a), with a the lattice con- 
stant. For a doping the YRZ model is described by a 
propagator, 

G[k,u:,x)=Y.'^^l\^-K-^lJ{^ + K)l (1) 

where = ±-Bfc , = y^^' + A^g, ^fc = (Cfc + 

^2)/2 and = ± , where gt{x) weights the 

coherent part.'* The energy dispersion = — 2t(cos k^a^ 
cos kya)—4:t' cos fc^a cos kya—2t" {cos 2kxa+cos 2kya)—iip 
includes hopping out to 3'^'^ nearest neighbour, while 
^1 = — 2i(cos k^a + cos kyO) is the first nearest neighbour 
term, which determines the placement of the pseudogap 
off the Fermi surface, coinciding with the AFBZ bound- 
ary. These energy dispersions contain doping dependent 
coefficients: t{x) = gt{x)to + 'igs{x)Jx/^: t'i^) = gt{x)tQ, 
and t"{x) = gt{x)tQ, where gt{x) = 2x/(l-|-x) and 
9s{x) = 4/(1 + x)'^ are the Gutzwiller factors. The dis- 
persion here uses fip as an effective chemical potential 
or Fermi level at T = 0, determined by the Luttinger 
sum rule. Values of other parameters in the dispersion 
were taken from Ref.^ to be: t' /to = -0.3, t" /to = 0.2, 
J /to = 1/3, and x — 0.338; while the optimal supercon- 
ducting gap Ag^ was chosen to give an effective optimal 
Tc around 90K for a ratio of 2A0^(r = 0)//cbTc = 6,i* for 
units of to = 0.1 eV, where A'^^{T) is the gap amplitude 
given in BCS theory at temperature, T. 

From the YRZ propagator of Eq. dT]), one can ex- 
tract the YRZ spectral function and see that there are 
four energy branches, given by the energies, ±£^^, where 
E"^ = \/ {E^y + Agj.. One can straightforwardly calcu- 
late the specific heat, C{T), from the entropy, S, given at 
temperature, T, by the standard formula summed over 
the four energy branches which reduces to 

S = -2kB if(^sc) HfiKc)] 

a=± k 

+f{-E^,)Hf{-E::)]}, (2) 

where / is the Fermi function and fee, the Boltzmann 
constant. The temperature dependence enters through 
both / and the temperature-dependent superconducting 
gap. The electronic specific heat gamma, denoted by 
7(r) = C(T)/T, is a constant, -/{T) = ^klN{0), in 
the non-interacting case, with N{Q), the electronic den- 
sity of states at the Fermi level. 7(T) is presented in 
Fig. [U for several doping values, normalized by the con- 
stant 7o = ^fc|iV(0,x — 0.16). These results are based 
on the generic phase diagram, illustrated by Fig. [H^a), 
which is slightly modified from that previously used by 
YRZ, with optimal doping now at x = 0.16. This is 
more relevant for comparison with the specific heat data 
of Loram et al^J^ Other phase diagrams based on the 
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FIG. 2: (Color Online) (a) The phase diagram via the mag- 
nitudes of the gaps A°^{x)/to = 0.24(1 - 82. 6(j: - 0.16)^ ) and 
Apg{x)/to — 0.6(1 — x/0.2). Also, the values of the critical 
angles, 6^^'^ and 6^^'^ vs x. (b) One quadrant in fc-space, 
identifying the critical angles where pseudogap exists. The el- 
lipse is a Luttinger pocket from the YRZ theory for x — 0.05. 



analysis of a large database have also appeared in the 
literatur o^^i^" , in which the pseudogap line ends near the 
upper edge of the superconducting dome, rather than 
a,t X = 0.2. Since, in this work, we are only interested 
in making a qualitative comparison with experiment, we 
have made no attempt to alter other basic parameters of 
the pseudogap state introduced in the original paper of 
YRZi in order to improve quantitative agreement with 
experimental data. 

Comparison of Fig. [T](a) with the experimental results 
of Loram et al^s^!^ shows that our theoretical results cap- 
ture all essential qualitative features observed. First, as 
one proceeds towards the underdoped regime, there is a 
significant decrease in 7(T), as Tc is approached from 
above. This reflects an effective decrease in the DOS 
around the Fermi energy in the normal pseudogapped 
state. Second, the jump at Tc is greatly reduced with 
increased pseudogap. In Fig. [Hb), we summarize the 
doping dependence of the normalized jump, AC/j{Tc)Tc, 
as well as the normalized condensation energy, AU/ AUn 
at T = 0, derived from our entropy calculations. Re- 
call that the internal energy, U, is related to the spe- 
cific heat by dU/dT — C{T) and ultimately to the en- 
tropy. The condensation energy AU is defined as the 
difference between U in the superconducting state, and 
its value in the normal state at T = 0. This is given by 
AU — Jq" {Snormsd — Ssc)dT . AC/n is the Condensation en- 
ergy when the pseudogap is set to zero in both the normal 
and superconducting states. Both normalized quantities, 
the jump and the condensation energy, are seen to drop 
precipitously with decreasing x due to the increase in 
pseudogap, which is not part of any pure BCS formula- 
tion, in which both quantities would be constant for all 
doping. Our findings agree qualitatively with the data of 
Loram et al^n^. It is clear that the model of YRZ has 
captured an additional essential element of the physics 
of the underdoped cuprates not present in standard BCS 
models. Note that our values of AC/'j{Tc)Tc are larger 
than experiments indicate and that this is mainly due to 
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FIG. 3: (Color online) Plots of 7(r)/7o comparing the super- 
conducting state (red dashed) to the pseudogapped normal 
state (solid), with the areal displacement shaded: (a) YRZ 
model and (b) nodal liquid. The arc model fitted to YRZ for 
a; = 0.13 and compared to the nodal liquid: (c) pseudogapped 
normal state and (d) the superconducting state. 



our use of a large gap ratio of 6, on the order of that in- 
dicated by STM.^* However, a smaller value would give 
better quantitative agreement, indicating that the STM 
data may not reflect the bulk. 

Next, we make connection with the nodal liquid and 
the Fermi arc approaches. For both cases, we assume 
that the pseudogap is located on the usual large Fermi 
surface [shown as the solid black curve in Fig.l^^b)]. This 
corresponds to replacing the AFBZ energy ^ of the YRZ 
model by the energy ^fc. In this limit, the expression of 
Eq. ^ reduces to the standard BCS expression, with the 
square of the superconducting gap replaced by the sum 
of the square of superconducting and pseudogap. For the 
Fermi arc model we apply an additional constraint that 
the pseudogap is non-zero only in an arc located around 
the antinodal direction, such that 



A" 

pg 



(cos k^a - cos hyo), {6 < 6*^^^) 



AO 

— ^(cos k^a 



cos kytt) 



> 



(3) 

and A^^*^ — 0, otherwise. The angle 9 — 
arctan [(tt — ky)/ (tt — kx)] is shown in Fig. El^b). For our 
purposes, 

^ARC 

is a single fitting parameter, adjusted to 
give 711(1" = 0) (n is for the normal pseudogapped state) 
equal to that of the YRZ model at the same doping. Ad- 
ditionally, we ignore any temperature dependence of the 
pseudogap magnitude below Tc as is suggested by Refs^^^ 



and"'^^ wherein they find that Raman data shows that 
there is little modification to the pseudogap below T^. 
We further ignore complications of the disappearance of 
the pseudogap at some, relatively high, temperature T* 
which is known to cause large specific heat anomalies at 
T* , which are not observed in experiment. In Fig. [2^a), 
the open red squares are the values of 9^^'^ obtained 
from the fit to the specific heat using the construction 
of Fig. lUb). In the YRZ model, the Luttinger pockets 
also form an ungapped region analogous to the Fermi 
arc model and it will be this region which is responsible 
for the specific heat. Consequently, for YRZ, we define 
9j^^ as the angle from (tt, tt) to the edge of the Luttinger 
pocket shown in Fig. [2l^b) and plot this in Fig. ^a,) as 
a function of doping in comparison with 0^^*^. The arc 
model fits consistently show 9^^^^ < 9j^^. This corre- 
sponds to the additional states (albeit with less quasipar- 
ticle weight^jA) which are located along the AFBZ bound- 
ary in the YRZ model, which are absent in the arc model. 
One might imagine unravelling the Luttinger pocket onto 
the Fermi surface of the arc model for a conceptual pic- 
ture of this comparison. Thus, we expect the arc model 
to capture much of the same features as the YRZ model 
for the specific heat. This contrasts to the nodal liquid 
case which has the gap over the entire Fermi surface and 
no fitting parameters. 

In the top two frames of Fig. [31 we compare supercon- 
ducting (red dashed curve) and normal pseudogap (solid 
black curve) results for 7(T) in the case x = 0.13. Frame 
(a) is for YRZ and frame (b) for the nodal liquid. The 
shaded areas illustrate the entropy difference between 
these two states. The entropy readjustment is less for 
the nodal liquid and consequently the specific heat jump 
at T = Tc is reduced. These differences arise because the 
density of states at the Fermi surface, N{0), is finite in 
the normal pseudogap state of the YRZ model while it 
is zero in the nodal liquid because the pseudogap exists 
over the entire Fermi surface, and thus, 7n(T' ^ 0) — * 
(n is for the normal pseudogapped state). A further com- 
parison of these two cases is presented in Fig. ^c) and 
[3l[d), where we have chosen to compare directly the two 
normal states and the two superconducting states, re- 
spectively. The nodal liquid and YRZ agree well at large 
T, but deviate significantly for T < 0. 025^0, with the 
nodal liquid curve going to zero. The open red squares 
are results for the arc model with 9^^^^ chosen to fit 
the value of 7^^^(T = 0). This results in a finite DOS 
at the Fermi level, for both the YRZ and arc models, 
which results in excellent agreement over all tempera- 
tures. Similar results for the superconducting state are 
presented in Fig.[3ljd). The agreement between YRZ and 
the arc model is excellent. Both show some slight devia- 
tions from the nodal liquid, but these deviations appear 
less important than in the normal pseudogapped state 
results of Fig. [3{c) where we find the low temperature 
differences to be quite striking. However, it should be 
noticed that, in Fig. [3]Jd), the nodal liquid result (dot- 
ted) does fall below YRZ for all T < Tc which results in 
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FIG. 4: (Color Online) (a) Condensation energy, AC/(0)/to, 
vs X for all three models, overlayed with Asc/^o and Apg/to 
from Fig. [2] (b) Plot of normalized jump vs x for all three 
models. Inset shows A7(rc)/7o. 



substantial loss of area under these curves, indicating a 
smaller condensation energy of the superconducting state 
[shown in Fig.[3]Ja)]. It is clear from this comparison that 
placing the pseudogap on the Fermi surface, rather than 
on the AFBZ, while at the same time cutting it off at 
^ARC rcpi-oduces well the YRZ results. For YRZ, there 
is no cut off, but rather the Luttinger contours move 
away from the gapped AFBZ boundary in the region of 
the nodes as is clear in Fig.[2l^b). 

Fig. m provides further comparison of both the nodal 
liquid and Fermi arc models with the YRZ results, across 
the entire doping range. Fig. Ufa) compares the conden- 
sation energy AJ7 vs doping, x. There is excellent agree- 
ment between YRZ and the arc model, while the nodal 
liquid is consistently lower. Although the nodal liquid 
still captures the large decrease in condensation energy 



that is caused by the opening of the pseudogap, it overes- 
timates the effect. The same remarks apply to Fig. Sfb) 
where AC/C„ is shown vs x. In the inset, we show the 
A7(Tc)/7o for completeness. 

In summary, we have found that the microscopic model 
of YRZ, based on the RVB spin liquid, which includes as 
its central essential element the formation of a pseudo- 
gap on the AFBZ in the underdoped region of the cuprate 
phase diagram, can account for all of the qualitative char- 
acteristics of the observed evolution of the specific heat 
as a function of doping. We have also found that an arc 
model with pseudogap formation on the Fermi surface 
itself, but limited to a region around the antinodal direc- 
tion resulting in an ungapped arc beyond 9c in the nodal 
region, can adequately simulate the results obtained in 
YRZ theory. The basic reason for this fortunate cir- 
cumstance is that, in YRZ, the Luttinger surfaces define 
pockets in the nodal direction around, but not directly 
on the antiferromagnetic Brillouin zone boundary, which 
keeps the electrons in that nodal region from fully sam- 
pling the pseudogap on the AFBZ. We have shown that 
this effect can be well approximated by the ungapped 
Fermi arcs of the arc model. On the other hand, the 
nodal liquid idea with pseudogap over the entire Fermi 
surface accentuates the pseudogap effect as compared to 
YRZ. Nevertheless, such a model still has merit; because 
of its great simplicity, it can give straightforward insight 
into the qualitative behaviour of the pseudogapped state. 
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